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(N ■ 

^ ■ Methods are presented allowing a realistic, experiment-based simulation of low- 

mass quark pairs coupling to a virtual photon. Differential cross-section reweighting 
factors allow to reproduce the resonant structure observed near quark pair produc- 
tion thresholds, and an algorithm is proposed that allows hadronization of the light 
quarks in the very low mass region, where existing theory-based models are ex- 
pected to break down. The corresponding routines can be embedded in two- and 
four-fermion generators for e^e~ machines. 



> 

m 
(N 



1 Introduction 

Several physics subjects of interest require a correct description of low mass quark- 
• antiquark systems. In e~^e~ qq processes, initial state radiation can bring back the 
mass of the outgoing quark pair near its production threshold, giving rise to topologies 
^ ! with an energetic photon and a collimated hadronic system. Such events can for example 
^ I be used to extract a measurement of the hadronic cross-section at the effective center- 
of-mass energy (and the ratio R of hadronic to pointlike muonic cross-sections) [1], or 
represent a background to searches for physics beyond the standard model, such as radia- 
tive decays of intermediate vector bosons. 

In neutral-current e~^e~ — > qqff processes, low-mass quark pairs give again access to 
hadronized virtual photons. Such events are interesting in their own right (they have 
been studied by the LEP experiments at LEPl [2], and the corresponding cross-section 
is measured by the Delphi experiment at LEP2 [3]), but can also represent a background 
to new particle searches, of which low mass Higgs boson searches are an example. 

Current two- and four-fermion generators are not adapted to these kinds of studies. A first 
point is that the string model [4], upon which the PYTHIA program [5] relies to simulate the 
hadronization of parton systems, describes the data well at high energy, but is expected 
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to break down below about 2 GeV. Secondly, computations of amplitudes generally make 
use of state-of-the-art knowledge of electroweak corrections and include a first order QCD 
term, but ignore non-perturbative corrections to vertices involving a light qq pair. In other 
words, diquark systems are generated with distorted mass distributions in the low mass 
region, and are hadronized unrealistically. 

Low mass quark- antiquark pairs appearing in e^e~ annihilations most often involve a 
coupling to a virtual photon. This suggests using experimental data on e"'"e~ — * hadrons a 
low energy, namely the R-ratio and measurements of exclusive final states cross-sections, 
to solve the problem. The mass spectrum and hadronization aspects are described in turn 
in the next two sections. Although it is not the main subject of this report, some results 
for two- and four-fermion processes are summarized in the third section. The code location 
is given in the last section. 



2 Quark-pair production near thresholds 

The task of implementing non-perturbative QCD corrections to qq differential cross- 
sections was first attempted in the framework of the four-fermion generator FERMISV [6,7]. 
The present method is similar and improves on a few shortcomings : the integration of 
the qq are spectra are improved using proper Monte-Carlo techniques ; most importantly, 
recent data in the charm resonance region are used in order to obtain a good description 
of this mass range. 

2.1 Reweighting of the 7*qq vertex 

With QCD turned off, the cross-section for quark pair production through a virtual photon 
can be written (in a symbolic way): 

Kadrons = Kn + ^^^dd + ^^4 + ^4 + ^^bb 

where the are just charge, colour and threshold factors, already taken into account in 
generators, and is the total parton level R-ratio. We now consider QCD effects, and 
write similarly: 

df^hadrons dff X R 

— da^^ X [Ru + Rd + Rs -\- Rc -\- Rb) 

= dal^R^ + dal^Rd + da%Rs + da^^R, + cia°j^i?6, 

with Rq = Rq/Rg .The above lines express that for a given quark ffavour q, a QCD- 
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improved description can be implemented by multiplying the parton level matrix elements 
by a factor Rg, representing the contribution of quark flavour q to the experimental R- 
ratio, with charge, colour and threshold factors divided outQ Cross-sections of processes 
containing a qq pair will thus be computed as 

da = ^g(mqq) X |Xp X dPh, 

with I Alp the parton level matrix elements, and Rg evaluated at the qq mass. 

The parametrizaton of R includes the contributions from all J^*" = 1 resonances up 
to the T(11020), with parameters taken from [8]. Intuitively (and as in [7]), the p's and 
tu's are attributed to the u- and d-quarks {ie, entered in the definitions of Ru and Rd); 
the 0's, the ip's and the T's are attributed to s-, c- and b-quarks (included in Rs,Rc,Rb) 
respectively. 

The continuum is parametrized by simple functions satisfying partonic boundary condi- 
tions, with parameters adjusted in the threshold regions so that the available data [9] 
are correctly reproduced. At present, one such function is defined for u- and d-quarks, 
one for s-quarks, and one for c-quarks. Since the currently available data do not allow a 
detailed representation of the Upsilon resonance region, the continuum from b-quarks is 
represented by a step function at the B-meson threshold. 

Finally, the continuum parametrizations include an overall second order perturbative QCD 
correction factor. The running of the strong coupling constant is computed at second order, 
and Aqcd is set so as to reproduce the world-averaged as{m%) [8]. The implementation 
used is taken from [5]. 

The resulting parametrization of the R-ratio is illustrated in Figure 1. 
2.2 Monte- Carlo sampling 

The R function being very accidented, special attention has to be paid to its Monte-Carlo 
sampling. Exploration of all narrow peaks is performed using iterative dividing of the 
integration range. In a first step, the full range is divided in a number of sub-intervals, 
the borders of which are set on the (a priori known) peaks of the resonances. The interval 
contributing most to the integral is then divided into two equal sub-intervals. This process 
is iterated and fine granularity is obtained by requiring a large number of intervals. After 
completion of the initialization phase, drawing intervals at random (in proportion to their 
contribution to the total integral) then ensures correct population of all the peaks. The 
program implementing this algorithm and used for our purpose, VESKO, was first described 
in [10]. It was modified to be able to take into account the positions of a priori known 
peaks in the distributions, and to handle several distributions simultaneously. 

^ Of course, the splitting of R into contributions from different flavours is only a computational 
convenience, and the whole procedure makes sense only after summing all contributions. 
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Fig. 1. Parametrization of the R-ratio, and comparison to various data sets. The upper figure 

ilhistratcs the range betweeen and 10 GeV. The data sets are from ND (below 1.5 GeV), BES 
(from 2 to 5 GeV), Crystal Ball (from 5 to 7.4 GeV) and MD-1 (between 7.4 and 10 GeV) [9]. 
The middle figure details the (p region, and the lower one shows the ^ resonance region. 

2.3 Caveats 

The above procedure is valid only if the final state quark pair is genuinely coupled to 
a virtual photon. Therefore, this method should not be applied to the e^e~ e~^e~qq 
process, unless a kinematic region is selected where the multi-peripheral (mostly 7-7) 
diagrams have negligible contribution and interference with the other diagrams. 

Other complications appear for final states with four quarks. In principle, one R fac- 
tor should be present for every pair of identical quarks, but there are cases where care 
should be taken. For example, in the uddu final state and if all quarks have same colour, 
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one should evaluate whether charged current diagrams (ie, WW production) or neutral 
current diagrams (Z/7* Z/7*) dominate the total amplitude, depending on momentum 

configurations. Similarly, in final states with four identical quarks (of same colour), the 
"pairing" which dominates the total amplitude has to be determined, and the R factors 
should be evaluated at the masses of the connected quark pairs. 

On the contrary, the e'^e~ — > qq process is unambiguous. 



3 Hadronization 



Having decided on the qq mass spectrum in the resonance region, it remains to be de- 
termined how every individual event should behave in the detector. As a function of the 
mass of the quark pair, it has to be decided whether we it can be treated as a resonance 
of known decay properties, and if not, if the string or cluster fragmentation algorithms 
available in PYTHIA will work. In the first case, the quark pair is replaced by the corre- 
sponding resonance in the PYTHIA event record, and it is decayed by a subsequent call 
to the hadronization routine. In the second case, it is assumed that above a mass of 2 
GeV, string hadronization is convenient, except in the cc and bb threshold regions where 
cluster hadronization correctly takes over. 

Since the lowest mass region is dominated by the p, we are left with a mass window 
between ~ 0.9 and 2 GeV, open to u-, d-, and s-quarks, where new input is needed. This 
input can be taken from experiment, since the first hadronic final states in e"'"e~ collisions 
have been extensively studied in the seventies and eighties. Data on the following exclusive 
non- resonant cross-sections [11] were compiled for this purpose: 

• e+e^ — >• Stt (7r"'"7r~7r''), 

• e+e- 47r (27r+27r-, 7r+7r-27r°), 

• e+e- Stt (27r+27r-7r°), 

• e+e- 6% (27r+27r-27r°,37r+37r-), 

• c+c- ^ 2K {K+K-,KsKl), 

• e+e- ^ KsK+TT- + CC , 

• e"'"e~ — > K'^K~TT'^7r~ . 

The KlK~^tt~ final state has never been measured due to experimental difficulties related 
to Kl detection, but was assumed to be equal to the KsK~^t:~ cross-section. 

These cross-sections define the "branching fractions" of the non-resonant 7*, which are 
displayed in figure 2. In practice, when a qq pair is found to be neither a resonance, nor 
heavy enough to be administrated by string fragmentation, one of the above final states 
is drawn at random according to the measured cross-sections at the current qq mass. 
The momenta of the outgoing hadrons are then generated uniformly in the available 
phase space using the CERN routine GENBOD [12] (it has currently not been attempted to 
reproduce the inner kinematic structure of these multiparticle final states. For example, it 
is ignored that a significant part of the 27r+27r^ final state is genuinely pp decays) . Finally, 
the quark pair is replaced by its daughter hadrons in the PYTHIA event record. 
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Fig. 2. The non-resonant {p and (p subtracted) branching fractions of the 7* as a function of its 
mass, from measurements of e"'"e~ — > exclusive final states [11]. 
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4 Results for specific processes 



The present algorithm has been implemented in the 2-fermion generator KK2f [13]. Run- 
ning this generator in the configuration of current e^e~ colliders allows to estimate the 
radiative hadronic cross-sections expected at these machines. Table 1 summarizes results 
for the existing $- and B-factories, and for a hypothetical Charm factory such as the Cleo-c 
program at CESR. Given the expected luminosities at these machines, (9(10^ — 10^) events 
are expected in regions where the knowledge of the hadronic cross-section is still limited[^. 
Figure 3 illustrates the differential hadronic cross-sections expected when operating at D- 
and B-meson threshold. 

In the context of four-fermion processes, events with a light quark-pair will show up as 
thin, collimated jet (most often a tt^tt' pair from a p decay) recoiling against an energetic 
lepton- or jet-pair. Table 2 summarizes some results for the /x"'"/i~qq and qqqq processes. 
Assuming 500 fb~^ recorded per experiment, the number of expected events with one 

^ These numbers are given with initial state radiation computed to third order (see [13]), and 
final state radiation switched off. A discussion on the level of accuracy of this approximation is 
beyond the scope of this work, but is expected to be 1-2%. 
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Table 1 

Expected e'^e~ qq event rates (in nanobarn) for various colliders. Cross-sections are given for 
three ranges of the effective centcr-of-mass energy (for the hypothetical Charm factory operating 
on the ^'(36'), the last number corresponds to the range 2-3 GeV). No selections are applied. 

Colhder ^ (GeV) 0-0.95 GeV 1.05-2 GeV 2-5 GeV 



DA$NE 1.02 42.38 

Charm factory 3.77 1.31 
KEKB,PEP-II 10.58 0.18 



0.77 
0.09 



0.95* 
0.25 



Fig. 3. Differential radiative cross-sections expected at the ^'(3S) resonance (left) and at the 
T(4S) resonance (right). No selections are applied. 
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quark pair of mass smaller than 2 GeV is 0(10^ — 10^), depending on the process under 
consideration. This confirms that such events can represent a substantial background to 
new particle searches, and need to be accounted for in the simulation of standard model 
processes. 

Table 2 

Expected e+e" /x"'"//~qq and qqqq event rates (in picobarn) at LEP2 energies. Quark-pair 
masses are restricted above 2 GeV (case 1), and released down to threshold (case 2). 

Final state ^/s (GeV) Case 1 Case 2 



qq 



qqqq 



200 
200 



0.33 
8.72 



0.36 
9.05 
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5 Program availability 

The routines implementing the methods described above are obtainable from [14] , together 
with installation notes and commented example applications. 



6 Conclusions 

A complete treatment of low-mass quark pair production and hadronization is described, 
which is relevant whenever the pair is produced via a virtual photon. This allows to 
release traditional cuts on quark pair invariant masses, and to study low mass hadronic 
systems specifically. The large statistics available at existing e~^e~ machines will allow 
to compare these predictions with experiment, with applications both in two- and four- 
fermion processes. 
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